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Abstract. 

Gamma-ray and electron-positron pair production will figure prominently in laser-plasma 
experiments with next generation lasers. Using a Monte Carlo approach we show that 
straggling effects arising from the finite recoil an electron experiences when it emits a high 
energy photon, increase the number of pairs produced on further interaction with the laser 
fields. 



1. Introduction 

Next generation lasers such as ELI [[Tl |2l and VULCAN [|3l IH are expected to achieve 
intensities of 10^^ W cm^^ or greater, opening up rich opportunities to study QED processes 
in the strong field regime[l5l [6l |71. One of these processes is pair production in counter- 
propagating beams, which is a significant effect at lO^^W cm~^|[8l, and, at lO^'^Wcm^^ is 
sufficiently strong that a cascade may develop in which each electron produces a pair that in 
turn produces further pairs, to initiate an avalanche. 

In intense laser beams, the dominant route to pair production is a two-stage process. 
Firstly, an electron with Lorentz factor y interacts with the laser fields {E,B) to produce a 
gamma-ray with high energy hv. Secondly, the gamma-ray interacts with the same fields as 
it propagates to produce an electron-positron pair |[8l|9l. The crucial parameters for the two- 
stage process are (i) rj « yl-E^ + vx Sl/Ecnli where ^cnt = 1.3 x 10'^ Vm"^ is the Schwinger 
field, and (ii) x = hv\E^ + (ck/k) x B\/2msC^Ecnt where hk is the gamma-ray momentum 
and Ej_ is the component of the electric field perpendicular to the direction of motion of 
the electron or photon, as appropriate. Strong field QED effects such as pair production 
become important when these parameters reach unity. The energy of the gamma-ray photon 
is typically hv « OAAriymc^, implying « 0.22i]^, and the cross-section for pair production 
is proportional to exp(-4/3;^f) for;^' ^ 1- Because of the exponential cut-off to the cross- 
section, the process is very sensitive to the value of rj reached by the electron, when the 
regime of strong QED is approached. This leads to an abrupt laser-intensity threshold for 

i This expression is generally accurate for y » 1, a precise definition is given in (|2]i 
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pair production as noted by [HI, S and [[TDl . where a small increase in t] can lead to a large 
increase in pair production. 

In lO, flU and [fTT]|. a semi-classical treatment was used in which electrons were subject 
to a continuous loss of energy through radiation of gamma-rays. In this continuous model, 
the radiation reaction force is approximated, to the lowest order in 1 /y, by the expression 



where g(r]) 6 [0, 1] accounts for the reduction of the total power radiated by the electron. 
F{ri,x) is the quantum synchrotron emissivity, as defined in Appendix A of [|9l. However, 
in reality, gamma-ray emission is a random quantum process [fT2l [T3l . As 77 approaches 
unity, the energy of the emitted gamma-ray becomes a substantial fraction, typically OAArj, 
of the kinetic energy of the electron before emission, and the resulting electron trajectory 
starts to fluctuate substantially away from the semi-classical average. This enhances pair 
production in two ways. Firstly, if an accelerating electron propagates by chance over an 
unusually large distance before emitting a gamma-ray, it may reach a large Lorentz factor y 
and emit a gamma-ray with an unusually high energy. This, in turn, has a much increased 
probability of producing a pair, as compared to the photons emitted by an electron moving 
on a classical trajectory in the same field. Such events cumulate, and we show in 13.11 that 
this straggling eff"ect W2\ enhances pair production. Secondly, the deviations in the trajectory 
caused by straggling lead the electron to sample the laser fields at locations other than those 
on the classical path. In [9], it was shown that, under the influence of continuous radiation 
losses, electrons naturally migrate towards points in the laser field at which they emit few 
high-energy gamma-rays. In the case of counter-propagating circularly polarised waves, 
these are the nodes in the electric field. Once an electron settles at an = node, its 
Lorentz factor drops and pair production ceases. In contrast, we show in 13.21 that straggling 
causes electrons to migrate more slowly towards these nodes. Consequently, electrons 
undergoing discontinuous momentum changes retain a large Lorentz factor for longer than 
their semi-classical counterparts, and emit a larger number of high energy gamma-rays, which 
subsequently produce more pairs. 

In this paper, we extend the calculations of |[8l and ^ by modelling radiation loss as a 
discontinuous process that leads to straggling effects on the electron trajectory. This is done 
using a Monte-Carlo method that is described in section [21 In section |3l we compare the 
continuous with the more realistic discontinuous loss case, and show that the latter produces a 
larger number of pairs at laser intensities around 10^^ W cm"^. Pair production at 10^"^ W cm~^ 
is relatively little affected by straggling since, at this intensity, even semi-classical, continuous 
electron trajectories lead to gamma-rays that are well over the threshold for pair production. 

2. Description of the Monte Carlo algorithm 




(1) 



In 4- vector notation, the quantities 77 and x are defined as 

eh i_ ,.i eh^ i_ 
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for an electron and photon whose four-momenta are p"" and Hk'', respectively, with F^y the 
electromagnetic field tensor. In the following, we make use of the approach presented in 
that retains the weak, quasi- stationary field approximation. The probabilities for gamma-ray 
photon emission and the subsequent pair production are in this approximation functions of rj 
and X only. An equivalent system is then chosen, for the same given x and 77, in which the 
electron moves in a plane perpendicular to a uniform, static B field. This allows us to use the 
transition probabilities as formulated in [[T4ll . 

The radiation emitted by the particle due to its acceleration in the laser fields — here 
called synchrotron radiation — is assumed to be a random walk process [fT2l [T3]| . The 
probability of emission of a gamma-ray photon is governed by an optical depth Tg. At the 
start of the calculation, and immediately following the emission of a photon, the current 
optical depth Tg is set to zero and a randomly chosen 'final' optical depth Tg^ is assigned to the 
electron. The current optical depth Tg increases as the electron propagates until it reaches r^, 
at which point a photon is emitted. At each computational timestep the electron's position, 
momentum and current optical depth are updated according to the equations 

— = -eiE-\-vAB) , (3) 
at 

at ntsY 

dTg r''d'N^ ^, 



where 



dxdt 

d^N V3 rj FiTj,x) _ 

dxdt luTc y X 

is the differential rate of production of photons of parameter x by an electron of parameter 
rj and F{ri,x) is the quantum synchrotron emissivity, expressed in the weak, quasi-stationary 
field approximation lfT4ll . An expression for F(ri,x) is given in Appendix A of BH. The 
Compton time is tq = ^c/c, where Ac = hlm^c is the Compton wavelength, and Off is the fine 
structure constant. 

When Tg = Tg^ is reached, a photon is emitted, and its energy is randomly assigned as 
follows. First, the parameter is found from the relation 



Jo dydt^ 



^ ^ - dxdt 

ri/2 djv dy 
Jo dxdt^ 



Jo X J 



(7) 



where ^ is a uniformly distributed random number in [0,1], and the right hand side is a 
monotonic, increasing function of;^f*^*\ which is tabulated in both the;^^^ and 77 directions. The 
variables and 77, in equation (|7]), depend on the electromagnetic field at the point of photon 
emission. Then, assuming that the photon is emitted parallel to the electron momentum, its 
energy hv,is, determined by the relation 

^^^2m^_ (8) 
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which follows from the definitions of 7] and x and the assumption that the electron is 
relativistic, 7 » 1 . For computational efficiency, we ignore the rare photons emitted in the low 
energy part of the spectrum, < x^^in' "^^^^^^min chosen such that the neglected photons 
carry only 10"^ of the total energy and are incapable of producing pairs. The computational 
procedure was successfully tested by its ability to reproduce the correct synchrotron spectrum. 
This is shown in Fig{TJ for a large number (10^) of emitted photons in a regime where 77 = 0.1. 
A x^'^^ parameter, which represents the photon energy, is assigned to each photon at the 
emission point, according to the sampling relation (|7]). The occurrence (in percent of the total 
number of the 10^ emission events) of emissions in each energy interval gives a distribution 
along the;^^^^ axis, that scales with the quantum synchrotron emissivity F{ri,x'"^^)- 
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Figure 1. Comparison between the distribution inx^^ resulting from the samphng relation ^ 
and the quantum synchrotron emissivity F(ri,x^^). 



After the electron has emitted a photon, which is assumed to occur instantaneously, it 
continues its trajectory beginning at the same point but with a different energy y®mc^ and 
momentum To ffnd p®, we again use the approximation that the photon is emitted parallel 
to the pre-emission electron momentum, so that, by conservation of momentum. 



P 



hv\ 
1 \p. 

cpj 



(9) 



In reality, the photon may be emitted over a cone of opening angle approximately I/7, but we 
neglect this small angular spread here. Equation ^ implies that a small amount of energy is 
extracted from the laser fields during photon emission, which we also neglect. This neglect 
is unimportant in the present test-particle model in which the laser fields are prescribed, since 
the fraction of the electron energy transferred to the laser field is small ~ 1/(77*^^), as shown 
in [Appendix A[ However this point may need more careful treatment in a more complete 



dynamic treatment of the development of an electron-positron cascade and its inclusion in a 
particle-in-cell code lfT5l[T6l . 
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3. Straggling effects 

As in [9J, we adopt the configuration of two counter-streaming, circularly polarized laser 
beams, which can qualitatively be considered equivalent to the interaction of an incident and 
reflected beam from a solid |[8l. We restrict the treatment to pair creation by the two-step 
process described in section [2l Pairs can also be produced directly by an electron interacting 
with the laser field, without the need to produce an intermediate real photon. This process 
is analogous to trident pair production by an electron in the Coulomb field of a nucleus. 
However, although this direct, single-step process may be important at lower laser intensity, 
it is dominated by the indirect two-step process at laser intensities above 10^^ Wcm~^, and we 
do not include it in the results presented here. 

In the two-step process, the emitted gamma-rays travel through the plasma and may 
produce pairs according to the cross-section which depends on;^. In many ways, the process 
is comparable with that of gamma-ray production by electrons. The photon trajectory could be 
followed through the electromagnetic fields and the optical depth integrated as for electrons. 
However, in order to make a direct comparison with the rate computed assuming continuous 
losses, we here follow the much simpler approximate approach that was used in [|9l. The 
probability of pair production is calculated for a gamma-ray with x = Af^*^ that propagates 
through electromagnetic fields that are assumed to equal those at the point where the gamma- 
ray was emitted. We follow [|9l and make the assumption that conversion into pairs can occur 
until the photon has propagated a distance of one laser wavelength. Here the approximation is 
that the beams have a transverse structure with cylindrical geometry, whose radius is one 
wavelength, representative of beams which are tightly focus sed at the diffraction limit to 
achieve maximum laser intensity. The real photons only produce pairs while they propagate 
through the region of large laser field. In contrast to electrons, whose trajectory is determined 
by the electromagnetic fields, the photons are undeflected. This means that all along its path 
through the laser fields, the value of the parameter remains unchanged atx = X^^- Further 
details can be found in section 3.3 of [l9l. The differential pair production rate via real photons 

d^^ 

is the product of the differential rate for real photon emission, (n^X)^ effective 

dxdt 

conversion probability into pairs, (l - exp(- (^p))). where (vp^ is an effective optical depth 
to pair conversion. The rate of pair production via real photons per electron is then addressed 
with the same rate equation as in [|9l 



'real 



dt 



^^d^^'^'^Hl-^^P^-^pW • (10) 



3.1. Cumulative effects of the straggling at a B = node 

In this section, we select a particular trajectory at a S = node (zero magnetic field) of an 
infinite standing wave generated by two counter-propagating, circularly polarized laser beams 
with the same sense of rotation of the fields. Throughout, distances and times are scaled with 
the laser wavelength and laser period, respectively. This configuration has been extensively 
analysed in |[8l |9l in the case of continuous radiation losses where it was shown that the 



Monte Carlo calculations of pair production in high-power laser-plasma interactions 



6 



electron settles into a circular orbit, which is favourable for pair production. Here we show 
how discontinuous trajectories increase pair production. 

The two beams counter- propagate along the z-axis, each with an intensity I24 x 
10^"* W cm'^^, with I24 = 0.3. The electron is initially located at a B = node in the standing 
wave formed by the counter-propagating beams. Its initial Lorentz factor is y = 10, with its 
momentum oriented along the positive j-axis. The electron motion is tracked over seven laser 
periods during which it remains at the magnetic node. In Fig. Wi^), the discontinuous Monte 
Carlo calculations lead to the straggling of the Lorentz factor for the particle (green curve). 
This is compared with the Lorentz factor of a particle that undergoes continuous radiation 
loss, computed as in @ (red curve) for the same initial conditions. A feature of the straggling 
in this regime is a ~ 30% fluctuation about the Lorentz factor derived from the continuous 
description. The particle motion, as shown in Fig. [2tb) is characterized by excursions to 
higher and lower 77 values than those reached by an electron undergoing continuous losses. 
The discrepancies in the particle y position, between the two descriptions, is also shown in 
Fig.|2tb). Finally, the comparison of Fig. [31 and Fig.[2tb), shows that the high 77 excursions are 
very effective at producing pairs. Moreover, the accumulation of pairs, produced by the high 
rj excursions, leads to an enhancement of the overall pair production, in this example by about 
40% after only seven laser periods. At this time, the average number of pairs produces by an 
electron reaches unity, which can be considered an approximate threshold for the initiation of 
an electron-positron avalanche. 
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Figure 2. Case of an electron at a B = node of two very long, counter-propagating, 
circularly polarized pulses. The trajectories found with continuous and discontinuous loss- 
descriptions are compared in terms of (a) Lorentz factor, (b) 77-parameter & spatial coordinate 

y- 



Monte Carlo calculations of pair production in high-power laser-plasma interactions 



7 




5 10 15 20 25 30 35 40 45 
time (laser phase) 



Figure 3. Number of pairs produced against time, by the continuous-loss trajectory (red curve) 
and the discontinuous one (green curve). 

3.2. Convergence properties of discontinuous trajectories, ojf-B = nodes 

In [|9l, the authors showed that the electron trajectory at a B = node is unstable. Small 
displacements from the node grow, causing the electron to migrate towards the stable E = 
nodes. When the electron reaches an E = node, its oscillatory momentum damps towards 
zero and pair production ceases. Here, we show that in the more realistic Monte Carlo 
description, the migration progresses less smoothly, and, on average, the electrons tend to 
take longer to reach to the E = node. This increases the number of pairs produced, since, 
on average, electrons emit high energy photons for a longer period of time. 

We consider the same configuration and initial conditions used in section 13. 1[ with 
I24. = 0.3, except that the initial electron position is now z = 0.1 instead of z = 0, which is the 
location of the B = node. Then the motion of a set of lO'* electrons is followed until the time 
^max = 20. Because of random fluctuations, each of the 10"^ trajectories is diff"erent, resulting 
in a diff"erent number of pairs finally produced. In Fig. Ufa), the z-position for two sample 
trajectories is tracked, and compared with the z-position of an equivalent electron starting 
from the same point but subject to a continuous radiation-reaction force. Two diff"ering kinds 
of behaviour are highlighted. In one case (the blue curve), an electron trajectory described 
by discontinuous radiation losses reaches the stable E = node earlier than the trajectory 
described by continuous radiation losses (red curve). In another discontinuous-loss case 
(green curve), the electron takes longer to converge on the stable E = node. The majority 
of electrons belong to this latter class. This is demonstrated quantitatively in Fig. Ilfb), 
where the red curve shows the evolution of the square of the distance to the E = node 
at z = n/2, averaged over 10^ discontinuous trajectories, i.e., ^(z - n/2)^y For comparison, 
the green curve shows the case of continuous radiation loss. The longer migration time, 
together with the excursions to large 77 discussed above, combine to enhance pair production 
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even further. The overall increase in pair production is shown in Fig. Sfc) which presents 
the distribution function for pair production compiled from 10"^ trajectories. The A^reai axis is 
divided into 100 equally spaced intervals in the logarithmic range logio(A'reai) € [-5, 1]. All 
the electrons are identically initialized at z = 0.1, with I24 = 0.3 (red bars). We find that the 
spectrum is substantially shifted upwards with respect to the continuous value (green line), 
logiQ(A/^,eai) - -1.01, obtained with continuous radiation loss. The figure shows that 77.5% 
of the particles produce more pairs than if their trajectory was computed in a continuous 
manner. A small number of the sampled trajectories come close to producing one pair per 
electron (logio(A^reai) - 0). 

For the same configuration and initial conditions, we now repeat the calculation at a 
lower laser intensity, I24 = 0.1. The results are presented in Fig. [2a), where the A^reai 
axis is again divided into 100 equally spaced intervals, this time in the logarithmic range 
logiQ(Meai) 6 [-6,0]. We find that the distribution is again shifted upwards with respect to 
the value logjo(A^reai) - -4.19 obtained with the continuous radiation loss. Electron energy 
straggling and longer migration times prove to have a significant impact on the number of 
pairs produced. For the top 1% of the trajectories, one order of magnitude more pairs are 
produced than are predicted using continuous radiation losses. 

In contrast, at higher laser intensities the spread of the distribution is quite narrow, and 
is centred close to the continuous value. This is illustrated in Fig. Ob), where we present 
results for I24 = 1. In this figure, the A^reai axis is divided into 100 equally spaced intervals 
in the logarithmic range logjQ(A/^reai) 6 [-2,4]. The location of the peak in the spectrum 
agrees well with the number of pairs computed with a continuous semi-classical description 
logiQ(A^reai) - 0.55. At high laser intensity, straggling has less effect, because the process is 
above threshold and less dependent on statistically rare excursions to large 77. 

In order to obtain a more straightforward figure of merit, in terms of produced pairs, we 
compute an effective value for the number of pairs produced per electron whose motion is 
followed: 

iel 

where / stands for the complete set of equally spaced intervals that discretize the A^reai axis. 
The index / refers to one of these intervals, in which the probability is P{i). This procedure 
amounts to summing over histograms such as those in Fig. ID^c), [2a) & [S^b), for the total 
expected number of pairs produced per electron starting from an initial position z = 0. 1 . 

For a range of laser intensities. Fig. [6ta) shows a comparison between the number of 
pairs produced per electron undergoing discontinuous losses (computed using a Monte Carlo 
approach), with the number found for continuous losses. For low laser energy, I^a = 0.1, the 
number of pairs (A^reai) associated with discontinuous-loss trajectories is about 5 times larger 
than the number associated with continuous-loss trajectories. This is due to the importance 
of stragglers in this regime, and has a strong dependence on the 77-parameter. Because of the 
large range in A^^eai (eight orders of magnitude) this effect is more easily seen in Fig [6tb), 
which plots the ratio between these two numbers, showing that (Nreai) /A'reai shrinks as the 
laser intensity increases, until it becomes very close to unity, at I24 = 1. The strength of the 
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Figure 4. The discontinuous and continuous descriptions are compared at a laser intensity 
/24 = 0.3, for initialization at z = 0.1, away from the B = node (located at z = 0). In (a) z is 
shown as a function of time, in (b) the square of the distance to the stable E = node (located 
at z = n/2), and in (c) the statistical distribution of the number of pairs A^reai produced by a 
trajectory. 



Straggling process is found to be very sensitive to the laser intensity for about I24 < 0.4. 
4. Conclusion 

In this paper, we report on calculations of the number of electron positron pairs produced as 
an electron moves in the field of a laser with intensity in the range range 10^^ - 10^"^ W cm~^. 
In this regime, the dominant route to pair production involves the emission of gamma-ray 
photons via synchrotron radiation, and their further conversion into electron-positron pairs 
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Figure 5. The statistical distribution (red bars) of a set of 10 electrons identically initialized 
at z = 0.1 (close to the B = node at z = 0) for (a) I24 =0.1, and (b) I24 - 1. The number of 
pairs produced by an electron undergoing continuous losses is shown for comparison (green 
lines). 




(a) (b) 

Figure 6. The continuous- and discontinuous-loss descriptions are compared in terms of the 
number of pairs produced, with respect to the laser intensity parameter I24, for two counter- 
streaming, circularly polarized pulses, the particle being initialized at z = 0.1, close to the 
B = node at z = 0. The curves are fits with rational functions. 
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E 13 • The new aspect of our treatment is that we take account of the stochastic nature of 
gamma-ray production, which leads to discontinuous quantum jumps in the electron energy; 
using a Monte-Carlo algorithm we follow the evolution of a large number of stochastic 
trajectories that result from each individual set of initial conditions. 

Pair production has been shown to be important in configurations with counter- 
propagating laser beams, that can be generalized to apply to laser-solid interactions |[8l lU. 
However, these authors used a continuous approximation for the radiation reaction force, and 
examined, therefore, one classical trajectory for each set of initial conditions. In order to 
make a direct comparison with our calculations, we investigate a specific, highly idealized 
configuration of the laser fields that was also used as a test-case by |[8l |9|l: that of two long, 
counterpropagating, circularly polarized laser pulses. 

Discontinuous, stochastic jumps in the electron energy give rise to a spread in energy 
as a function of time known as straggling. Our main result is that straggling increases the 
number of pairs produced for given initial conditions, and for a given laser intensity. This 
behaviour is more pronounced for laser intensities less than 0.4 x lO^'^Wcm"^ where the 
transition probabilities are strongly dependent on rj and x\ in this range we observe up to a 
fivefold increase in the number of pairs produced. Straggling efi'ects are less important at 
higher intensities, where most gamma-rays are above threshold for pair production and the 
number of electron-positron pairs produced is large. 

Two efi'ects contribute to the increase in pair production due to straggling. Firstly, the 
value of the electron's 77-parameter, which controls the energy of the emitted gamma-ray, 
varies considerably about its mean. Occasional excursions to large 77 result in the emission 
of higher energy gamma-rays with a correspondingly larger probability of producing pairs. 
Secondly, straggling increases, on average, the time taken for electrons to congregate in a 
region where few pairs can be created (a node with £^ = in our calculations), and this extends 
the time during which high-energy photons are emitted and pairs produced. We conclude that 
straggling lowers the laser intensity required for significant pair production. 

In realistic situations, when large numbers of pairs are produced, their mutual interactions 
as well as their back-reaction on the accelerating fields must be taken into account. In an 
astrophysical setting, [|T6| has addressed this problem by combining a Monte-Carlo algorithm 
with a PIC-code. We envisage that the Monte Carlo algorithm we present will serve as a 
basis for a coupled Monte Carlo-PIC algorithm aimed at modelling the initiation of prolific 
pair-plasma creation in intense laser fields. 
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Appendix A. Estimation for the energy transferred to the laser field by a single particle 

From the momentum change expression at the photon emission time p^*^ = p-lik, we estimate 
the error, in the Monte Carlo algorithm, related to the energy transferred to the laser field by 
a single particle. The expression for the relative error reads 

Ay 1 I «^ 2 I 

— = —\y'-y-\- hv/(meC ) 

7 7 

In the limit where the 1 /y and 1 /y^^ parameters are small, it can be approximated as 



Ay ^ 1 1 / 1 
y 2y\y*^ 



7 



1 



) 



which is of the order 1 /(yy^^). 



